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We study the effect of nonmagnetic impurities in the superconducting state of the t — J modeL 
The Bogoliubov de Gennes equation is derived using a slave-boson mean-field approximation, 
and we solve it numerically at T = 0. The d^2_y2-wave order parameter is suppressed near an 
impurity and its spatial variation induces the extended s-wave component. We also find the local 
charge density wave (CDW)-like feature due to the Friedel oscillations. This leads to the splitting 
of the expected zero-energy Andreev bound state. Therefore the system is unlikely introducing a 
localized state of other symmetry such as superconductivity with broken time reversal symmetry 
or antiferromagnetism. 
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The symmetry of the superconducting (SC) state in 
high- Tc superconductors (HTSC) has been a subject of 
intensive study since it is an important clua to clar- 
ify the mechanism of their superconductivity.EP Now it 
is estabhshed that the SC state has a predominantly 
da;2_j,2-wave character with a possible mixture of an s- 
wave component due to the orthorhombic lattice distor- 
tion in some systems. In d-wave superconductors the 
effect of nonmagnetic impurities is quite different from 
that of conventional (s-wave) superconductors because 
of the phase structure of the pair wavefunction and the 
resulting presence of nodes in the excitation gap. This 
problem has been studied intensively bTt|ixi|ami_authors, 
mostly based on weak coupling theories.nBcl'BB' In this 
letter we examine the effect of nonmagnetic impurities 
in the SC state of the t-J model on a square lattice. 
This model df;acribes the low-energy electronic proper- 
ties of HTSCu'EI' where strong correlation effects for the 
electrons are known to be important. Mean-field (MF) 
theories based on a slave-boson method nredjct a super- 
conducting state with a d-a-prf-symmetryBll^P They rnaif 
also explain the magneticU'tJ as well as the transports 
properties of HTSC if the gauge fields representing the 
fluctuations around the MF solutions are taken into ac- 
count. Therefore the effect of nonmagnetic impurities 
in thfi-SC states of the t-J model is of particular inter- 
est.0B 

We treat the t-J model with the Hamiltonian including 
that of the nonmagnetic impurities 
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where and I denote nearest- neighbor bonds and the 
impurity sites, respectively, and Cia is an electron opera- 
tor within the Hilbert space excluding double occupancy. 
Here Si is the spin-1/2 operator at a site i and J(> 0) 



is the antiferromagnetic superexchange interaction. The 
impurity potential Vq is taken to be Vb ^ J- We use 
the slave-boson method to enforce the condition of no 
double occupancy by introducing spinous {fia', fermion) 
and holons (6^; boson) {cia — blfia)- We decouple this 
Hamiltonian by a mean- field approximation (MFA). In 
the following we consider only the case of T = 0, so that 
holons are Bose condensed. Then the mean-field Hamil- 
tonian is written in terms of spinous only. 
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where S and A are the doping rate and the chemical po- 
tential, respectively, and fi = ii, ±y. The summations 
for i and j are taken over all sites. Here the singlet RVB 
order parameter (OP), Ay, and the bond OP, Xij^ 
defined as 



Ay = ihlfj] - fi].fjl)/2, Xl. 
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respectively, for each nearest neighbor bond. We con- 
sider the solution where Aij may be complex, while Xij 
is taken to be real (i.e., Xij = Xji)- We do not consider 
the possible antiferrromagnetic correlation here. 

Without impurities Hmfa can be easily diagonalized 
by the Bogoliubov transformation assuming spatially 
uniform x ^nd A. In the presence of inhomogeneities, 
however, we have to take their spatial variations into 
account. They are most conveniently described by the 
Bogoliubov de Gennes (BdG) equation.113' From eq.(2) 
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it is obtained as 
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Here {un{i),Vn{i)) is the wave function at a site i, and 
e„ is the corresponding energy eigenvalue. The OP's Xij 
and Aij, and the doping rate S (we take A as an input 
parameter) are expressed in terms of Un{i) and Vn {i), 
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where /(e) is the Fermi distribution function and N is 
the total number of lattice sites. 

We numerically solve the self-consistency equations (5) 
and (6) . For simplicity we assume that Xij has a spatially 
uniform value (obtained by MFA without impurity) ex- 
cept for the bonds connected to the impurity sites where 
Xij = 0. First we choose A, and assume some initial 
values of A^- and 5. Inserting these into the BdG equa- 
tion we diagonalize the resulting matrix using House- 
holder method (the size of the matrix is 2M^ x^M^(= 
2N X 2N), if the system size is M x M(= N))\l3 Then 
we recalculate Ay and 5 according to eq.(6). This pro- 
cedure is iterated until the convergence is reached. From 
Ay we define a d^2_y2- and an extended s- wave OP 
component on the site i, 
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where x and y correspond to one lattice constant in x- 
and y-direction, respectively. We have checked that the 
values obtained for the system without impurities agree 
well with those in the usual mean-field calculation, when 
M > 16. In this case only Ad is finite and it has a 
spatially uniform value. 

Now we turn to the results for the system with an im- 
purity. We consider the case where the concentration 
of impurities is low so that the states around the im- 
purities can be treated independently. We fix t/ J — 3 
throughout in the following, and the doping rates are 
chosen to be 5 = 0.20. (We have studied the system 
with 0.05 < 5 < 0.20. The results are qualitatively the 
same for all cases.) In Fig.l the spatial variation of A^ is 
shown, which is suppressed near the impurity, and the ef- 
fect is strongest along the diagonals of the square lattice. 
This is due to the interference effects for momenta close 
to the gap nodes (sign change of the pair wave function). 
The extended s-wave component is induced in the region 
where A^ is not uniform (Fig. 2), and A,, at sites rotated 
90 degree around the impurity have the same magnitude 
but the opposite sign. It vanishes for sites located along 
the diagonals passing through the impurity site. {Ad and 



As are always real in the solutions we found.) 

We analyze the above results qualitatively by using 
the GL theory (though GL theory is not quantitatively 
valid at T = 0). The GL free energy icir the system with 
an impurity is written generically asH 
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where the last line represents the effect of the impurity 
located at r = (g^ > 0, gs > 0), and the crystal axis 
directions are denoted as x and y. We note that the cou- 
phng term like (A^A* +c.c.) should not arise due to the 
symmetry. Due to the gd term A^ is suppressed and its 
gradient becomes finite over the range of the coherence 
length. Then As is induced through the mixed gradient 
(-R')-term. In the K-ievm. the gradients in x and y direc- 
tions have opposite sign following the A^-symmetry, and 
the induced As must refiect this property. Therefore Ag 
should change sign under 90-degree rotation around the 
impurity. When the continuous change of the phases of 
Arf and As were allowed. As could change sign under 
90-degree rotation with keeping its amplitude finite (i.e., 
(pds = ±7r/2 along the diagonals). If the relative phase 
between A^; and As {(t>ds) is neither or tt, the state has 
{d + is)- or (tL|-pJs)-symmetry and breaks time reversal 
symmetry T.ESESP In this case the gap nodes disappear 
in the vicinity of the impurity and the system could gain 
more condensation energy. In our calculation no solution 
of this type with A^ and As in a complex combination 
appeared for the doping range 0.05 < S < 0.20. In the 
following we would like to discuss a possible reason for 
this result. 

It has been discussed previously based on the T-matrix 
approximation that non-magnetic impurities iii,a d-wave 
superconductor create a subgap bound stateja^ similar 
to Shiba's bound state in a ccuiventional superconductor 
around a magnetic impurity .Eir In the T-matrix formu- 
lation it is found that a single bound state occurs at 
the Fermi level (zero-energy) irLthe unitary limit, i.e. 
for a strong impurity potential.^ This is understood 
also ^Jiithin the BdG formulation as an Andreev bound 
state.Ed' Because a strong scattering center destroys the 
gap in the near vicinity of the impurity, quasiparticles 
can be trapped in this potential well and is subject to 
Andreev reflection at the walls of this well. Their en- 
ergies are determined by the phase of the gap function 
in the momentum directions which the particle-hole tra- 
jectories connect via scattering at the impurity .E3 For 
the d-wave superconductor there are two possible phase 
differences, and tt. The former leads to a state at the 
bulk gap value while the later generates a zero-energy 
bound state similar to the one at the , 
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face of a dj.2_j^2-wave superconductor .E^''^@SS*i311' 
The presence of the latter bound states should yield a 
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finite contribution of local density of states (LDOS) at 

the Fermi energy. In this case it could be energetically 
favorable to twist the OP phase in order to shift the 
zero-energy bound state away from the Fermi energy. 
Such a state would then correspond to the state with 
locally broken time reversal symmetry around the im- 
purity as mentioned above. This effect is, of course, a 
matter of competition between the energy gain of the 
quasiparticles and the energy expense due to phase gra- 
dients (yielding supercurrents) introduced by the twist. 

In order to examine this view we consider the LDOS 
which is denoted as N{ri,u) on the site i, 

N{ri,oj) = ^ [\un{i)\^5{uj - e„) + \vn{i)\^5{uj + e„)] . 

n 

(9) 

In Fig. 3 we show the result for N{ri,Lu) for a nearest 
neighbor site of the impurity {6 = 0.20). At an en- 
ergy much lower than the gap value (A^ ~ 0.12) we 
find bound states. The bound state energy is different 
from the naively expected value, i.e., zero, leading to a 
small gap between the state below and above the Fermi 
energy. The LDOS is different for these two states, since 
the doped system is not particle-hole symmetric. 

The splitting of the expected zero-energy bound state 
into levels slightly above and below zero occurs obviously 
without the T-violation, since our solution has only real 
superconducting OP's. This splitting may have various 
reasons. The validity of arguments based on the Andrccv 
approximation might not be guaranteed here. Another 
important aspect is, however, connected with charge den- 
sity oscillations as shown in Fig. 4. They originate from 
the Friedel oscillations due to the presence of the im- 
purity. The dominant Q- vector lies along the [110]- 
dircction which is for 5 = 0.2 very close to Q = (tt, tt) and 
leads to a nearly commensurate charge density stagger- 
ing. In the vicinity of the impurity these lead to a charge 
density wave (CDW) like environment with an effective 
"doubling of the unit cell" so that the local opening of a 
gap is expected. We would like to emphasize, however, 
that this CDW feature is not the result of an instabil- 
ity triggered by the quasiparticle density of states at the 
Fermi level, but is driven by the presence of the impurity. 
The fact that there is no quasiparticle peak in the LDOS 
at zero energy leads to the conclusion that a transition to 
a state with local T-violation or antiferromagnetic order 
is very unlikely to occur. 

In summary we have studied the effects of nonmag- 
netic impurities in the superconducting state of the t-J 
model using the Bogoliubov de Genncs equation derived 
via a slave-boson mean-field approximation. Near the 
impurity the d-wave OP is suppressed, and the extended 
s-wave component is induced as expected also from the 
GL description. The CDW-like feature due to the Friedel 
oscillations is most likely responsible for the double peak 
structure foimd in the LDOS at energies close to zero. 
We conclude that the violation of time reversal symmetry 
is suppressed here due to the removal of a quasiparticle 
bound state at zero energy. This is in contrast to the 
situation at the surface or at interfaces between d-wave 



superconductors. 
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Fig. 1 The spatial variation of in the system with 
an impurity. Here S = 0.20, t/ J = 3 and the system size 
is = 16 X 16. 

Fig. 2 The spatial variation of A^ in the system with 
an impurity. Here 5 = 0.20, t/ J = 3 and the system size 
is AT = 16 X 16. 
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Fig. 3 LDOS at a nearest neighbor site of the impu- 
rity with S = 0.20, tjj = 3 and A/' = 18 X 18. Here we 
have introduced finite width F = 0.008J to each state. 

Fig. 4 The electron density around the impurity. Here 
(5 = 0.20, tjJ = 3 and the system size is = 18 x 18. 
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